Why is the nuclear symmetry energy so uncertain at supra-sat urat ion densities? 
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Within the interacting Fermi gas model for isospin asymmetric nuclear matter, effects of the 
in-medium three-body interaction and the two-body short-range tensor force due to the p meson 
exchange as well as the short-range nucleon correlation on the high-density behavior of the nuclear 
symmetry energy are demonstrated respectively in a transparent way. Possible physics origins of 
the extremely uncertain nuclear symmetry energy at supra-saturation densities are discussed. 
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I. INTRODUCTION how to treat quantum many-body problems at high den- 

sities, the various techniques used in different many-body 
theories may be among the possible origins of the very 
The density dependence of nuclear symmetry energy uncertain E sym (p) at supra-saturation densities. Never- 



E sym (p) is currently a key issue in both nuclear physics 
and astrophysics, see e.g., Refs.[l-ll]. Despite much 
theoretical and experimental effort, our current knowl- 
edge about the E sym (p) is still rather poor especially 
at supra-saturation densities. Experimentally, some con- 
straints on the E sym (p) at sub-saturation densities have 
been obtained recently from analyzing nuclear reaction 
data, see, e.g., Refs. [12-14]. At supra-saturation densi- 
ties, however, the situation is much less clear because of 
the very limited data available and the few model anal- 
ysis carried out so far although some indications of a 
super-soft E sym (p) at high densities have been obtained 
from analyzing the ir + /ir~ ratio in relativistic heavy- 
ion collisions [15]. Theoretically, essentially all available 
many-body theories using various interactions have been 
used in calculating the E sym (p) . Unfortunately, the pre- 
dictions at supra-saturation densities arc very diverse, for 
a recent review, see, e.g., Ref. [6]. Assuming all models 
arc equally physical and noting that there is no first prin- 
ciple guiding its high-density limit, it is fair to state that 
the E sym (p) at supra-saturation densities is currently still 
completely undetermined. So, why is the E sym (p) so un- 
certain at supra-saturation densities? This is obviously 
an important question that should be addressed timely 
especially since several dedicated experiments have now 
been planned to investigate the high density behavior 
of the E sym (p) at CSR/China [16], GSI/Germany [17], 
MSU/USA [18] and RIKEN/ Japan [19]. Identifying the 
causes for the uncertain high-density E sym (p) may help 
experimentalists to decide what experiments to do and 
what observables to measure. While we can not fully an- 
swer this question, we identify several important factors 
and demonstrate their effects on the high-density behav- 
ior of the E 3ym (p) using probably the simplest many- 
body theory available, namely the interacting Fermi gas 
model for isospin asymmetric nuclear matter, e.g., Ref. 
[20]. There are many long-standing physical issues on 



theless, it is still very useful to examine effects of some 
common ingredients used in most many-body theories, 
such as the three-body and tensor forces, in the sim- 
plest model possible. While the interacting Fermi gas 
model can not be expected to describe all properties of 
infinite nuclear matter and finite nuclei as accurately as 
those more advanced microscopic many-body theories, 
it does give an analytical expression for the E syrn (p) in 
terms of the isospin-dependent strong nucleon-nucleon 
(NN) interaction in a physically very transparent way 
[21]. Most importantly, the key underlying physics re- 
sponsible for the uncertain E sym (p) at supra-saturation 
densities can be clearly revealed. In particular, effects of 
the spin-isospin dependent effective three-body force, the 
density dependence of the in-medium short-range ten- 
sor forces and the short-range nucleon correlation can be 
demonstrated clearly. The results are expected to be use- 
ful for not only understanding predictions of the various 
many-body theories but also ultimately determining the 
E S ym(p) at supra-saturation densities. 



II. SYMMETRY ENERGY WITHIN THE 
INTERACTING FERMI GAS MODEL 

According to the well-known Lane potential [22], the 
single-nucleon potential U n / p can be well approximated 
by 



U n / p (p, k) rj U (p, k) ± U sym (p, k)5 



(1) 



where the Uo(p,k) and U sym (p,k) are, respectively, the 
isoscalar and isovector (symmetry) nucleon potentials. 
Within the interacting Fermi gas model for isospin asym- 
metric nuclear matter [20], the nuclear symmetry energy 
can be explicitly expressed as (for detailed derivation of 
this formula, please see Ref. [21] and references therein). 



'Corresponding author, Bao-An_Li@Tamu-Commcrce.edu 



IP ( n \ — J?kin , rppotl 



J sym \ 

1 dt 

6 dk 



ppot2 
sym 



(2) 



\k F 



I dUp 
6 dk 



-jU sym (p, k F ), 



2 



where t(k) = h 2 k 2 /2m is the kinetic energy, m is the 
average nucleon effective mass and k F = (3tt 2 p/2) 1 / 3 
(Thomas- Fermi Approximation [23] ) is the nucleon Fermi 
momentum in symmetric nuclear matter at density p. 
We notice here that the Eq.(2) is identical to the one 
derived earlier by Brueckner, Dabrowski and Haensel 
[24, 25] using K-matriccs within the Brueckner theory. 
From Eq.(2), it is seen that the symmetry energy E sym (p) 
is only dependent on the single-particle kinetic and po- 
tential energies at the Fermi momentum kp. The first 
part E^y^ — g^(^§-)' P^ is the trivial kinetic contri- 
bution due to the different Fermi momenta of neutrons 
and protons; E^"^ = \k F -k F is due to the mo- 

mentum dependence of the isoscalar potential and also 
the fact that neutrons and protons have different Fermi 
momenta; while the E^°^ = \U sym (p,k F ) is due to the 
explicit isospin dependence of the nuclear strong interac- 
tion. The Uo(po,k) at normal nuclear density po is rela- 
tively well determined from the nucleon optical potential 
obtained from the Dirac phcnomenological model anal- 
ysis of nucleon- nucleus scattering data [26]. Moreover, 
interesting information about the U${p,k) at abnormal 
densities in a broad momentum range has been obtained 
from transport model analysis of nuclear collective flow 
in heavy-ion reactions [27]. For the momentum depen- 
dent part of the isoscalar potential Uo(k,p), we use here 
the well-known GBD (Gale-Bertsch-Das Gupta) param- 
eterization [28] 



Ugbd(p, k) = 



-75p/p 



1 + (k/(Ak F )Y 



(3) 



where A 



1.5. The El^l is then given by 



= 75p/ Pc 



V(3A 2 ) 
(1 + (1M 2 )) 2 



(4) 



The GBD potential describes reasonably well the 
nucleon-nucleus optical potential and has been widely 
used in transport model simulations of heavy-ion reac- 
tions [29]. Similar to the Ef y n ml the always in- 
creases with density. However, the U sym (p,k) is very 
poorly known especially at high densities/momenta [6]. 
To reveal the fundamental physics responsible for the un- 
certain high-density behavior of the E sym (p), we denote 
uto = u' as the n-p interaction in momentum-space in 
the isosinglet (T=0) channel, while Uti = u nn = u pp = 
u np is the nuclear strong interaction in the isotriplct 
(T=l) channel. In the latter, the charge independence 
of strong interaction has been assumed. Then the single- 
nucleon mean-field potentials are [30, 31] 
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Therefore, 

Uo{p, k) = -([/„ + Up) = -(3ut! + uto), (6) 
U sym (p, k) = w^{U n - Up) = - uto)- 

Thus, the U sym (p, k) measures the explicit isospin depen- 
dence of the nuclear strong interaction, namely, if the n-p 
interactions were the same in the isosinglet and isotriplct 
channels, then the U sym (p,k) would vanish. Currently, 
the calculation of U sym (p,k) is rather model dependent 
[31, 32]. In fact, its value can be either positive or neg- 
ative at high densities/momenta [33-35], leading to the 
dramatically different predictions on the high-density be- 
havior of the E sym (p). In coordinate space, in terms of 
the two-body NN interactions Vro( r ij) an d Vri(rij) for 
the isosinglet T=0 and isotriplet T=l channels, respec- 
tively, we have [20] 



E sym = \u sym (p,k F ) = -{V T \ 
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where 



Vto = tjP 



VTo{nj)d s rij, 
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is determined by the competition 

between the Vti and Vto, and thus by the isospin de- 
pendence of the in-medium NN interactions. The latter 
is largely unknown and is actually a major thrust of re- 
search at various radioactive beam facilities around the 
world. In fact, how the nuclear medium may modify the 
bare NN interaction and properties of hadrons has long 
been one of the most critical issues in nuclear physics. For 
instance, it is very difficult to obtain the empirical sat- 
uration properties of symmetric nuclear matter by using 
the bare NN interactions within non-relativistic many- 
body theories. Various in-medium effects, such as the 
many-body force [36, 37], the relativistic effect relating 
the in-medium interactions to free space NN scattering 
[38, 39] , or the in- medium tensor force due to both the 7r 
meson and p meson exchanges [40-42] , have to be consid- 
ered to obtain a reasonably good description of satura- 
tion properties of symmetric nuclear matter. How these 
effects manifest themselves in dense neutron-rich matter 
may affect the high-density behavior of the nuclear sym- 
metry energy. 



III. EFFECTS OF THE SPIN-ISOSPIN 
DEPENDENT THREE-BODY FORCE 



In this section, we examine effects of the spin-isospin 
dependent three-body force on the symmetry energy at 
supra-saturation densities. As it has been shown repeat- 
edly in the literature, see, e.g., Refs. [36, 43-45], a zero- 
range three-body force can be reduced to an effective 
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two-body force 



V d = t Q (l + x P tT )p a S(r), 



(9) 



where to, a and Xo are parameters and P a is the spin- 
exchange operator. Depending on how many and which 
properties of finite nuclei and nuclear matter are used in 
reducing the three-body force, various values have been 
used for both the xq and to parameters. Moreover, the 
density-dependence p a is often used to mimic additional 
in-medium effects such as the many-body force. For in- 
stance, in relativistic approaches, the dressing of the in- 
medium spinors will also introduce a density-dependence 
to the interaction, which can lead to similar effects as 
the many-body force. The parameter a can thus take 
on any value between and 1. The parameter xq con- 
trols the relative contributions to the three-body force 
from the isosinglet and isotriplet NN interactions. The 
three-body term has been included in many effective in- 
teractions, such as the Skyrme and Gogny interactions, 
in both phcnomenological, e.g., Rcfs. [46, 47] and micro- 
scopic, sec, e.g., Ref. [48], many-body theories for isospin 
asymmetric nuclear matter. Noticing that the potential 
energies due to the three-body force in the T=l and T=0 
channels are, respectively (for details, see Table II and 
the corresponding equation E ST in Ref. [45] ) , 
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one sees immediately that the terms containing xo cancel 
out in calculating the EOS of symmetric nuclear matter. 
In fact, among the normally 13 parameters used in the 
Hartree-Fock calculations with the Gogny force Xo is the 
only one having this special feature. To evaluate quanti- 
tatively the symmetry energy we use the Gogny central 
force [45] 



Vc(r) 



i=1.2 



(Wi + B,P a - H.Pr - M^Pr) 



2 / 2 



(11) 

where P T is the spin exchange operator and the values 
of the parameters W, B, H, M and p are taken di- 
rectly from Ref. [45]. It is necessary to stress here that 
since we are aiming at a qualitative understanding of why 
the symmetry energy is so uncertain at supra-saturation 
densities, the model parameters are not re-tuned self- 
consistcntly to reproduce any existing constrain on the 
symmetry energy at and/or below the normal nuclear 
matter density po- The resulting E?°™ is 



(i- 



-2x )^p a+1 . 



(12) 



Shown in Fig.l are the symmetry energy functions ob- 
tained with different values for the xo and a parameters. 
It is clearly seen that for a given value of the a parameter, 
the Xq controls the high density behavior of the symme- 
try energy. By varying the xo one can easily cover the 
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FIG. 1: (Color online) The symmetry energy obtained with 
different spin and density dependences in the three-body 
force. 



whole range of symmetry energy calculated within the 
Hartree-Fock approach using over one hundred Skyrme 
and Gogny forces [47] without changing anything in the 
EOS of symmetric nuclear matter. In particular, with 
xo = 1 and a = 1/3 as in the original Gogny force 
[45], then only the (S=l, T=0) spin-isospin n-p inter- 
action contributes to the three-body force and the sym- 
metry energy, the E sym (p) drops quickly to zero above 
certain supra-saturation densities. We regard this kind of 
E S ym(p) as being super-soft. Overall, depending on the 
value of the parameter xo , the symmetry energy can be 
cither stiff that keeps increasing with density or becomes 
super-soft above certain supra-saturation densities. 

In the 2003 survey by J. R. Stone et al. [47] of 
87 Skyrme interactions, the xq {x^ in the notation of 
Ref. [47]) ranges between —1.56 to 1.92. As a special 
example, we notice that B. A. Brown used xo between 
0.03 to 0.9 with the SKX [49]. The values of x used 
in Fig.l are thus within the uncertainty range of xq 
found in the literature. Moreover, we notice that the 
result here is consistent with the earlier finding that the 
EOS of pure neutron matter is very sensitive to the xq 
parameter [49, 50]. To our best knowledge, it is cur- 
rently not clear how to fix the xo alone experimentally. 
Nevertheless, it is worth noting that the cross section 
of pn charge exchange reactions and the symmetry 
potential U sym extracted from the isospin dependence 
of the nucleon optical potentials arc all directly related 
to the spin-isospin dependent nuclear interaction. They 
may thus be used to constrain the value of xq. Moreover, 
the density dependence of the nuclear symmetry energy 
itself can be used to constrain the xo should it be 
determined experimentally in the future. However, 
as we shall discuss in the next section, the in-medium 
tensor force may have similar effects on these observables. 
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IV. EFFECTS OF THE IN-MEDIUM 
SHORT-RANGE TENSOR FORCE AND 
NUCLEON CORRELATION 

Studies based on microscopic many-body theories indi- 
cate consistently that the symmetry energy is dominated 
by the isosinglet (S=1,T=0) channel [30, 51]. While mea- 
sured properties of deuterons indicate unambiguously 
that a tensor force is at work in the (S=1,T=0) chan- 
nel. It is also well known that the n and p meson ex- 
changes contribute to the intermediate-range attractive 
and the short-range repulsive tensor forces, respectively, 
according to [40] 



Vr(r) 



and 



-in 



-Tl-T 2 (-5l 2 )[ 
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The 



Su = 6(5 • r) 2 /r 2 - 2S 2 = 4S 2 P 2 {cos(6)) is the ten- 
sor operator. We notice that some effective interactions, 
such as the Paris force [52], only consider the tensor force 
due to the ir exchange and has a short-range cut-off. As 
shown in Fig. 3 of Ref. [53] , the short range behavior of the 
tensor forces used in popular effective interactions differ 
dramatically. Moreover, different short range cut-offs are 
normally introduced, e.g., in studying the single particle 
energy levels in rare isotopes [53, 54]. Unfortunately, the 
different behaviors of the tensor force at short distance 
being cut-out, thus not probed from studying the energy 
levels of single particles, affect dramatically the high den- 
sity behavior of the symmetry energy. Furthermore, in- 
medium properties of the p meson may affect the strength 
of the short-range tensor force. Because the symmetry 
energy is very sensitive to the competition between the 
isosinglet (Vro) and isotriplet (Vti) channels, the con- 
tribution of the in-medium short-range tensor force that 
exists only in the T=0 channel may thus affect signifi- 
cantly the high-density behavior of the E sym (p). To in- 
vestigate effects of the in-medium tensor force, we use 
the Brown-Rho Scaling (BRS) for the in-medium p me- 
son mass according to m*/m p = 1 — a BR ■ p/p [41]. 
Shown in Fig. 2 is the radial part of the total tensor force 
Vt = V£(r) + Vf(r) at p = 0,p o ,2p O7 and 3po, respec- 
tively, with the BRS parameter aBR = 0.2. As one ex- 
pects, the total tensor force becomes more repulsive in 
denser matter when the p meson mass is reduced accord- 
ing to the BRS. While the experimental evidence for the 
BRS is still not very clear, we use it here as an effective 
way of adjusting the in-medium strength of the tensor 
force. This is useful to mimic different ways of involving 
the tensor force in many-body calculations in the litera- 
ture. 

Noticing that the tensor force gives no contribution at 
the mean-field level to the potential energy with spheri- 




FIG. 2: (Color online) The radial part of the tensor force 
V T = V,f(r) + Vf{r) at p = 0, p , 2p , and3/9 with the BRS 
parameter aBR = 0.2. 
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FIG. 3: (Color online) The symmetry energy with different 
values of the BRS parameter a BR = 0,0.05,0.10,0.15,0.20 
using three values for the correlation parameter r\. 



cally symmetry tensor correlations [55], we assume here 
that the tensor force acting in the isosinglet n-p chan- 
nels in nuclear matter behaves the similar way as in 
deuterons, namely, the tensor operator S\ 2 is a constant 
of 2. Moreover, we introduce a two-step tensor correla- 
tion function, i.e. f(r) = 0, for r < r c and f(r) = 1, 
for r > r c where r c = ri(S / A-k p) 1 ' 3 is the "healing dis- 
tance" or short-range cut-off. Thus the tensor contri- 
bution to the isospin T=0 channel symmetry energy is 
Vr = J f(r)[V£(nj) + Vf{nj)]d a rij. The parameter rj 
is used to effectively vary the short range cut-off. The 
density dependence in the r c reflects the design that the 
size of nucleons shrinks as the density increase. With 
i] = 2, the two nucleons will always keep in touch on 
their surfaces while their sizes decrease with increasing 
density. 

Shown in Fig. 3 are the total symmetry energy func- 
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tions obtained using three typical values for the r\ param- 
eter. With rj = 10, the r c is so large that the tensor con- 
tribution to the EP°%% is completely cut off. The E sym (p) 
thus keeps increasing with density due to the and 
the EP^ terms. With smaller rj values, the high den- 
sity behavior of the E sym {p) is controlled by the BRS 
parameter ckbr- Similar to varying the xq parameter in 
the three-body force, by varying the cxbr parameter the 
tensor force can lead to E sym (p) from stiff to super-soft. 
A recent study indicates that a value of ocbr ~ 0.15 is 
required to reproduce the measured lifetime of 14 C [56]. 
With such an cxbr, the symmetry energy can easily be- 
come super-soft with ij between 0.1 and 1. Results for the 
three typical cases shown in Fig. 3 clearly indicate that 
the high-density behavior of the symmetry energy is sen- 
sitive to both the short-range in-medium tensor force and 
the NN correlation function. The short-range repulsion 
generated by the p-meson exchange plays the key role 
in determining the symmetry energy at supra-saturation 
densities. In fact, the relationship between the short- 
range repulsive tensor force in the isosinglet n-p channel 
and the appearance of the super-soft symmetry energy 
was first noticed by Pandharipande et al. within varia- 
tional many-body (VMB) theories [57, 58]. Ultimately, 
because of the dominance of the repulsive n-p interaction 
in isospin symmetric nuclear matter at high densities, it 
is possible that pure neutron matter is energetically fa- 
vored, leading to the negative symmetry energy at high 
densities. Of course, this can only occur if the p tensor 
contribution is sufficiently strong due to, for example, its 
reduced mass in dense medium. It was also pointed out 
that the fundamental reason for the completely differ- 
ent high-density behaviors of the E sym predicted by the 
VMB and the Relativistic Mean Field (RMF) models is 
the lack of the p tensor contribution to the energy in the 
RMF models [59]. 

To this end, it is interesting to note that the three- 
body force and the tensor force can affect similarly the 
high-density behavior of the symmetry energy. This is 
very similar to the situation in describing the saturation 
properties of symmetric nuclear matter. It has been 
shown that the saturation properties can be equally 
well described by either including the three-body force 
or the in-medium NN interactions based on the BRS 
[42, 48]. Nevertheless, as it was pointed out in Rcf. [44], 
since the three-body force is essentially a convolution of 
two-body forces, a consistent three-body force should 



also include a tensor component. We have only studied 
here separately effects of the three-body force and the 
two-body tensor force on the high density behavior of 
the symmetry energy. Effects including both the two and 
three-body tensor forces simultaneously will be investi- 
gated in a forthcoming work. In particular, by varying 
parameters controlling both the three-body force and 
the tensor force together we shall demonstrate how large 
the parameter space is in which the symmetry energy 
may become negative at supra-saturation densities. 

V. SUMMARY 

In summary, the high density behavior of the sym- 
metry energy E syrn has long been regarded as the most 
uncertain property of dense neutron-rich nuclear matter 
because even its trend is still controversial. Within the 
interacting Fermi gas model, the E sym can be expressed 
in terms of the isospin-dependent NN strong interaction. 
The high-density behavior of the E sym is determined by 
the competition between the in-medium isosinglet and 
isotriplet nucleon-nucleon interactions. Respective ef- 
fects of the in-medium three-body interaction and the 
short-range tensor force on the high-density behavior of 
the E sym are examined separately. It is found that the 
strength of the spin (isospin) dependence of the three- 
body force and the in-medium p meson mass in the short- 
range tensor force are the key parameters controlling the 
high-density behavior of the E sym . These findings are 
useful for understanding why the nuclear symmetry en- 
ergy is very uncertain at supra-saturation densities. 
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